A&A manuscript no. 

(will be inserted by hand later] 



Your thesaurus codes are: 

08 (08.02.2; 08.03.5; 08.13.2; 08.09.2; 13.21.5) 



ASTRONOMY 

AND 
ASTROPHYSICS 



The warm circumstellar envelope and wind of the G9 lib 
star HR 6902 



o 

o ■ 
o ■ 
<N : 

=3 : 

< 

^ ■ 
CNi : 



T. Kirsch, R. Baade, and D. Reimers 

Hamburger Sternwarte, Universitat Hamburg, Gojenbergsweg 112, D-21029 Hamburg, Germany 
Received ; accepted , 



> 

oo : 
i> 

m 

00 ' 

o ; 
o . 
o 

Oh; 

6 ■ 



Abstract. lUE observations of the eclipsing binary sys- 
tem HR 6902 obtained at various epochs spread over four 
years indicate the presence of warm circumstellar material 
enveloping the G9 lib primary. The spectra show Si iv and 
Civ absorption up to a distance of 3.3 giant radii (Rg). 
Line ratio diagnostics yields an electron temperature of 
^ 78 000 K which appears to be constant over the ob- 
served height range. 

Applying a least square fit absorption line analysis 
we derive column densities as a function of height. We 
find that the inner envelope (< 3 Rg) of the bright gi- 
ant is consistent with a hydrostatic density distribution. 
The derived line broadening velocity of ^ 70 kms^^ is 
sufficient to provide turbulent pressure support for the re- 
quired scale height. However, an improved agreement with 
observations over the whole height regime including the 
emission line region is obtained with an outflow model. 
We demonstrate that the common (3 power-law as well 
as a P oc p wind yield appropriate fit models. Adopting 
a continuous mass outflow we obtain a mass-loss rate of 
M = 0.8 - 3.4 X 10~^^ M0 yr~^ depending on the partic- 
ular wind model. 

The emission lines observed during total eclipse are 
attributed mostly to resonance scattering of B star pho- 
tons in the extended envelope of the giant. By means of a 
multi-dimensional line formation study we show that the 
global envelope properties are consistent with the wind 
models derived from the absorption line analysis. We ar- 
gue that future high resolution UV spectroscopy will re- 
solve the large-scale velocity structure of the circumstel- 
lar shell. As an illustration we present theoretical Si IV and 
C IV emission profiles showing model-dependent line shifts 
and asymmetries. 

Key words: binaries: eclipsing - binaries: spectroscopic - 
mass loss - Stars: individual: HR6902 - Ultraviolet: stars 



1. Introduction 

The spectroscopic binary system HR 6902 (G9 Hb -|- B8- 
9 V) belongs to the category of ( Aur stars, where a 
hot companion (generally a main sequence B star) moves 
through the extended envelope of an evolved late-type 
star. The specific binary geometry provides an exceptional 
opportunity to probe the spatial dependence of the phys- 
ical properties of the outer atmosphere and wind of the 
primary. Observations of ^ Aur systems and related ob- 
jects with the International Ultraviolet Explorer {lUE) 
and more recently with the Hubble Space Telescope (HST) 
take advantage of the ultraviolet spectral region. At wave- 
lengths below 3000 A the hot secondary dominates the 
observed flux allowing a separation of the spectral infor- 
mation. In the UV the observed spectrum is purely that 
of the companion with superimposed circumstellar and 
interstellar lines. Further information on the binary tech- 
nique can be found in sever al reviews (e.g ., Reimers 1989; 
iGuinan 1990l [Ahmad 1993[ [Harper 1996] [Baade 1998[) 



Ground based spectroscopic as well as UV and X-ray 
observations from space of stars in the cool half of the 
HR diagram have shown that the extended atmospheres 
of giants and supergiants undergo a dramatic physical 
change across a region in the HR diagram. To the right of 
this transition we observe cool and massive winds plus 
extended cool chromospheres, while to its left there is 
much hotter circumstellar matter. This transition is ac- 
companied by different spectral indicators defining dis- 
tinct dividing lines. The leftmost is the "X-ray dividing 
line". All stars to its left show X-ray emission, proba- 
bly arising from extended corona-like plasma (e.g., Haisch, 
Schmitt, & Rosso 1991, 1992). In the ultraviolet the same 
stars also show C iv and Si iv and often N V emission 
( Linsky fc Haisch 1979| ). Additional dividing lines are de- 
fined by asymmetries in the Ca ii and Mg ii emission lines 



(e.g., Stencel fc MuUan 1980a, 1980b). The strongest and 
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most reliable evidence for the existence of a cool, massive 
wind is the appearanc e of circumste llar absorption fea- 
tures in the Can line ( Reimers 1977 ). A theoretical dis- 
cussion of the dividing line concept based on hydromag- 
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netic waves can be found in Rosner (1993) and Rosner et 
al. (1995). 

However, this simple picture is complicated by the ex- 
istence of the so-called hybrid stars. Hartmann, Dupree 
& Raymond (1980) and Reimers (1982) found that some 
bright K giants and G supergiants show blue-shifted Mg ii 
(and partly Ca ii) absorption and emission lines like C iv 
and Si IV signifying temperatures of at least 10^ K. Hy- 
brid stars are also X-ray emitters ( Rcimcrs et al. 1996D . 
The outflow properties of these stars are still very un- 
certain and it is controversial whether hybrids show cold 
(~ 10^ K) or possibly warm 10^ K) winds. 

Previous lUE observations of HR 6902 suggest the ex- 
istence of warm circumstellar material that may repre- 
sent the transitional physics of hybrid star atmospheres. 
Reimers, Baade, & Schroder (1990a) and Reimers et al. 
(1990b) have presented a series of lUE high resolution 
spectra (i? 10*) taken a few days before and after 
eclipse. Si iv and C IV column densities indicate a height- 
independent electron temperature of « 70 000 K. As 
discussed by Reimers et al. (1990b) the density distribu- 
tion is compatible with a low mass-loss wind. 

In the present paper we report the results of a reanal- 
ysis of lUE spectra of HR 6902. Compared to the first 



Table 1. lUE observations of HR 6902. All exposures 
were taken through the large aperture 



study (Reimers et al. 1990a, 1990b) we utilize a substan- 



tially expanded observations database. The paper is or- 
ganized as follows: In Sect. § we present the observations 
and basic stellar parameters used in this paper. Sect. ^ 
summarizes the methodical principles of our absorption 
analysis technique and the application to the HR 6902 
spectra. In Sect. ^ we outline the multi-dimensional ra- 
diative transfer procedure used to infer global properties 
of the circumstellar shell. 

2. Observations and basic stellar data 

While nearly all observations have been conducted by our- 
selves in several lUE programs, the data used here are 
taken from the lUE Final Archive. We have compiled a 
sequence of large aperture spectra obtained between 1987 
and 1990. Full details of the NEWSIPS data reduction 
have been given by Garhart et al. (1997), and will not 
be repeated here. It is, however, worth reiterating that 
the new processing techniques generally greatly enhance 
the quality of lUE spectra. The increase in signal-to-noise 
improves the detection limit of faint absorption features 
significantly. Furthermore, the new noise handling allows 
a reliable error estimate for each point in the flux spec- 
trum. The archival lUE spectra exploited for this paper 
are listed in Table |^. 

The analysis presented by Reimers et al. (1990b) is 
based on lUE observations of the 1988 and 1989 eclipses, 
covering a range of impact parameters p (projected dis- 
tance between the stellar components) from ^ 1.1 — 2.0 gi- 
ant radii. Our new study also includes the high resolution 
spectra of the 1990 eclipse showing Si iv and C iv absorp- 



Date 


SWP 


Exp. time 


P 


Phase 




No. 


(min) 


(Rg) 




1987 










Jul. 11 


31325 


140 


5.26 


0.9227 


1988 










Aug. 25 


34131 


140 


1.15 


0.9887 


Aug. 25 


34135 


110 


1.10 


0.9898 


Sep. 03 


34177 


130 


1.19 


0.0122 


Sep. 04 


34182 


125 


1.33 


0.0147 


Sep. 05 


34187 


100 


1.48 


0.0175 


1989 










Jun. 30 


36589 


140 


9.12 


0.7926 


Sep. 26 


37191 


180 


1.68 


0.0209 


Sep. 27 


37195 


190 


1.84 


0.0235 


Sep. 28 


37200 


194 


2.00 


0.0261 


1990 










Sep. 20 


39668 


195 


3.32 


0.9534 


Sep. 20 


39669 


195 


3.29 


0.9538 


Sep. 23 


39699 


195 


2.81 


0.9612 


Sep. 23 


39700 


180 


2.79 


0.9616 


Sep. 26 


39713 


177 


2.28 


0.9694 


Sep. 28 


39721 


170 


1.98 


0.9742 


Oct. 08 


39795^ 


410 


0.84 


0.9997 



^ Low-dispersion spectrum 

tion up to p = 3.3 Rg. These additional impact parameters 
allow us to directly examine the wind acceleration region. 
One additional observation obtained 1987 yields an upper 
limit to the column densities at p = 5.3 Rg. 

A direct analysis of the circumstellar Si iv and C iv 
lines is not feasible owing to the considerable photo- 
spheric background absorption. We extract the circumstel- 
lar contribution by dividing all high resolution spectra by 
the out-of-eclipse exposure SWP 36589. This observation 
serves as suitable reference spectrum since the large im- 
pact parameter of p = 9.1 Rg implies a pure B star spec- 
trum. Indeed, the final wind model presented in Sect. |^ 
confirms that there should be no detectable circumstellar 
absorption. As an example we present in Fig. |l| the phase 
variation of the C IV absorption line. The observed line 
strength depends strongly on the relative position of the 
B star and decreases with increasing impact parameter. 
Though circumstellar Cil absorption at 1335 A is clearly 
detectable, the poor S/N ratio of the lUE spectra does 
not allow an adequate absorption line reconstruction. In 
the case of nearly saturated photospheric absorption lines 
the dividing procedure leads to spurious line profiles with 
artificial emission features. Unfortunately, the C iv and 
Si IV doublets are the only circumstellar absorption lines 
that can be used for a quantitative analysis. 

Great significance must be attached to the emission 
line spectrum visible during total eclipse of the B star. The 
1988 total eclipse observation (SWP 34143, (j) = 0. 



i/ = 0.9534 




((1 = 0.9538 




((1 = 0.9612 




= 0.9616 















1548 1550 1552 
Wavelength (A) 

Fig. 1. Series of lUE high resolution spectra of Civ 
UV mult. 1 and the corresponding positions of the B star 
behind the primary as seen from Earth. The spectra are 
smoothed with a Savitzky-Golay filter 



Table 2. Orbital elements and stellar parameters of HR 
6902 [gathered from Grifhn et al. (1995); the distance is 
obtained from the Hipparcos Main Catalog] 



Parameter 


Value 


Period (days) 


385.0 


Eccentricity 


0.311 


Inclination (deg) 


87 


Longitude of periastron (deg) 


146 


Semimajor axis (cm) 


2.94 X 10^^ 


Periastron passage (JD) 


2 447 306.33 


Conjunction (JD) 


2 447 788.18 


Spectral type (primary) 


G9 lib 


Spectral type (secondary) 


B8-9 V 


Rg/Rq 


33.0 


J?b/R© 


3.0 


Mg/Mq 


3.86 


Mb/Mq 


2.95 


Teff.G (K) 


4900 


reff,B (K) 


11600 


Mv,G (mag) 


-1.81 


Mv,B(mag) 


-0.03 


£;B-v(niag) 


-0.18 


Distance (pc) 


247 ± 40 



was part of the first HR 6902 program and has been re- 
ported by Reimers et al. (1990a). The deep 410™ exposure 
SWP 39795 (Fig. |) obtained on 1990 October 8 provides 
a significantly improved signal-to-noise ratio and will be 
subject of this paper. We note that the large difference 
in quality does not allow a comparative analysis in or- 
der to search for long-term variations. While in "normal" 
C Aur systems hundreds of scattering lines (e.g., Fell and 
Si ii) are seen, only a few are observed in the spectrum 
of HR 6902. A comparison with stars of similar spectral 
type suggests that the emission line spectrum cannot be 
due to the G star alone (Reimers et al. 1990a). Instead, 
the observed C iv and Si iv doublets are formed primarily 
by resonance scattering of B star photons in the circum- 
stellar envelope. 

Griffin et al. (1995) have redetermined the orbital and 
stellar parameters of HR 6902 by means of photometric 
observations during the eclipses of 1987 and 1989 (Ta- 
ble H). The relative orbit and the eclipse geometry are 
shown in Fig. ^. It should be noted that the stellar param- 
eters of giants belonging to C Aur systems can be much 
better determined than for single stars. A recent evolu- 
tionary analysis suggests that the bright giant primary 
of HR 6902 is already evolved into its blue loop phase 



(Schroder et al. 1997) 



3. Absorption line diagnostics 

A comparison of the total eclipse emission line spectrum 
with the out-of-eclipse observations shows immediately 
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Fig. 3. Relative orbit and binary geometry for all phases 
with detectable Si iv and C iv lines 



that the continuum flux is orders of magnitude larger than 
the scattered line flux. Hence, the circumstellar absorption 
profiles can be interpreted as pure absorption lines formed 
in the extended envelope of HR 6902. In the subsequent 



sections we present details of the analysis and interpreta- 
tion of the spectra of HR 6902. 

3.1. Basic equations 

The line formation in HR 6902 is dominated by velocity 
fields on different scales. To simplify the concept we in- 
troduce a line broadening velocity Wb containing both the 
thermal velocity uth and a stochastic component Wgto as- 
sociated with turbulence or other short-scale nonthermal 
motions. Assuming a Gaussian distribution of the stochas- 
tic velocity the spectroscopic line broadening parameter 
becomes = w^j^ -I- v'^^^. The macroscopic velocity field 
is characterized by a spherically symmetric expansion de- 
scribed by an a priori unknown velocity law Uwind('')- 

Throughout this paper we assume the line broadening 
velocity to be constant in the relevant parts of the enve- 
lope. It is convenient to define a dimensionless frequency 
variable x ^ [v — vq)//^.^^, where z/q is the line-center 
frequency and Ai^d — VQVy^/c the corresponding Doppler 
width. Macroscopic velocities are measured in the same 
units, i.e. u{r) = v{r)/vh. Regarding the B star as a point 
light-source of intensity Ib the emergent flux for a given 
impact parameter p can be written 



Fx{p) = Ib oxp [-Tx{p)] 



(1) 



The optical depth along a ray from the observer to the B 
star is given by 



n(r)$ 



-u{r) 



dz, 



(2) 



where denotes the dimensionless profile function. 

The geometry is defined by the distance from the primary 
r and the integration path z — — P^- The bound- 
aries Zmin and Zmax depend on the position of the B star 
and the adopted outer radius of the envelope (chosen as 
large that the results are unaffected by its actual value), 
respectively. The line opacity per absorber reads 



ninC 



(3) 



where / is the oscillator strength and Aq the rest wave- 
length of the transition under consideration. 

In the static case (i.e. without considering the expan- 
sion explicitly) Eq. ^ degenerates to a simple relation be- 
tween the optical depth and the column density N viz 



(4) 



Our analysis is carried out using a non-linear least 
square fit technique based on the Levenberg-Marquard 
method (e.g.. Press et al. I992D to derive the column den- 
sity N and the line broadening parameter Vsto- The 
NEWS IPS error estimate provides statistical errors of the 
deduced parameters and allows to assess the goodness- 
of-fit. As an additional constraint the best-fit parameters 
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Assuming the simplifying conditions of a low den- 
sity plasma the C iv and Si iv density ratio can be used 
to derive an ionization temperature Tion- We solve the 
ionization balance using the radiative-coUisional equilib- 



rium program CLOUDY (version 9003; Ferland 1997). 
The metal abundances are supposed to be solar {Aq = 
3.6 X 10-^,Asi = 3.6 X 10-5 jAndcrs fc Grevcssc 1989t ). 
As an approximation we identify the ionization temper- 
ature with the kinetic temperature knowing that the 
ionization state may freeze in, so that we have Tg ^ Tjon 
in the outer wind. It turns out that the electron temper- 
ature appears to be constant within the observed part of 
the envelope {p — 1.1 — 3.3 Rg) with a weighted mean of 
Te = 78 000 ± 2000 K. However, the ion ratio may be af- 
fected by different CNO abundances due to the first dredge 
up which can deplete the carbon abundance significantly 
(e.g., Lambert 1981 ). For example, a carbon abundance 
reduced by a factor of 4 would lead to a temperature of 
- 90 000 K. 



Fig. 4. Column densities of C iv and Si iv obtained by the 
quasi-static approach 

are forced to be equal for all lines of the particular mul- 
tiplet under consideration, i.e., we favor the strategy of 
resonance doublet modeling instead of analyzing individ- 
ual absorption lines. Our fitting procedure allows to vary 
the continuum flux level in order to avoid the difficulties 
associated with a continuum definition by eye. Finally, 
the theoretical profile has to be convoluted with the lUE 
line-spread function. The instrumental profile is assumed 
to be well represented by a Gaussian with FWHM in ac- 
cordance with the actual determination of the NEWSIPS 



resolution ( Garhart et al. 1997 ) 



3.2. Quasi-static approach 

As a starting point of our analysis we use a quasi-static 
approximation neglecting macroscopic velocity fields. The 
objective of this approach is to provide a provisional set of 
column densities and stochastic velocities. However, line 
broadening effects due to the wind outflow calls the static 



model into question. In Sect. 3.3 we will examine possible 



consequences for the derived parameters. 

In the spectrum SWP 31325 (impact parameter p = 
5.2 Rg) we cannot identify any signature of Civ or Siiv 
absorption lines. Instead, we use the error information of 
the lUE flux to determine upper limits of the column den- 
sities. The results of the static analysis are summarized in 
Table |. Fig. | shows the run of column densities as a 
function of the impact parameter. As already noticed by 
Reimers et al. (1990b), the egress observations indicate 
significant differences from epoch to epoch. The obtained 
stochastic velocities show a considerable scatter with a 
weighted mean of 68 ± 16 kms"^. 



3.3. Density models 

With known electron temperature we are able to trans- 
form the ion density into total hydrogen densities. In the 
next step we try to match the derived hydrogen column 
densities employing different types of density laws. Each 
density law implies a physical model defined by a specific 
set of parameters. Again we implement a least square fit 
procedure using the Levenberg-Marquard method. 

3.3.1. Hydrostatic model 

As a first attempt we use the hydrostatic equation for an 
isothermal plasma 



nu[r) = no exp — - 

r H 



(5) 



where rg defines the bottom of the envelope with the cor- 
responding hydrogen density no- H denotes the density 
scale height. This approach is motivated by the fact that 
at least the subsonic region will be well represented by the 
hydrostatic density law. In order to fit the empirical col- 
umn densities we adjust the density parameter no and the 
scale height H. Identifying ro with the giant radius the 
least square procedure yields no = 5.0 ± 0.92 x 10^ cm^'^ 
and H = 3.A± 0.14 x 10" cm. 

The pressure may be composed of two parts, the ther- 
mal gas component and an effective turbulent pressure 
leading to the total dynamical pressure 



( YT 



''turb 



(6) 



where is the mean molecular weight measured in units of 
the hydrogen mass. The velocity fturb denotes an isotropic 
turbulent motion with a spatial scale smaller than the 
pressure scale height. Assuming the ionization equilibrium 
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Table 3. Column densities, stochastic velocities, and temperatures derived from least square profile modeling. The 
parameters are obtained using two distinct wind models. At each phase the first line displays the results of the quasi- 
static approach. The second line shows the effect of the macroscopic velocity field (/3 power-law with Voo — 100 kms^^). 
Additional information can be found in the text 



Phase 




log iV [Si IV] 






Ic 


g N[C IV 


] 




i'sto(kms "*") 


r(io^K) 






logA/'p] 









9227 






< 12.95 








< 13.30 










— 

















9534 


13 


27 


(-1-0.11, 


-0 


15) 


13 


60 


(+0.12, - 





16) 


101.9 (+ 23.1) 


76.4 


(±5 


0) 


18 


08 


(-l_0.18, - 





11) 






13 


25 


(-H0.10, 


-0 


13) 


13 


59 


+0.10, - 





13) 


62.7 (± 29.4) 


76.4 


(±5 


3) 


18 


06 


(+0.15, - 





12) 





9538 


13 


24 


(+0.12, 


-0 


16) 


13 


56 


(+0.10, - 





13) 


64.0 (± 23.0) 


76.4 


(±4 


5) 


18 


04 


(-1-0.18, - 





12) 






13 


23 


(+0.11, 


-0 


14) 


13 


59 


(+0.10, - 





13) 


65.5 (± 31.7) 


77.3 


(+5 


6) 


18 


04 


(+0.17, - 





14) 





9612 


13 


31 


(+0.11, 


-0 


14) 


13 


90 


(+0.07, - 





09) 


103.3 (± 16.2) 


82.2 


(±3 


2) 


18 


16 


(+0.19, - 





12) 






13 


32 


(+0.10, 


-0 


13) 


13 


91 


(+0.06, - 





07) 


74.0 (+ 20.3) 


82.2 


(±3 


7) 


18 


17 


(+0.16, - 





16) 





9616 


13 


28 


(+0.10, 


-0 


13) 


13 


72 


(+0.08, - 





10) 


58.5 (± 12.5) 


78.9 


(±3 


2) 


18 


10 


(-1-0.15, - 





11) 






13 


30 


(+0.10, 


-0 


12) 


13 


77 


(+0.08, - 





09) 


48.8 (+ 22.3) 


79.4 


(±4 


1) 


18 


13 


(-1-0.16, - 





16) 





9694 


13 


61 


(+0.05, 


-0 


06) 


14 


07 


(+0.05, - 





06) 


64.2 (± 5.9) 


79.4 


(±1 


9) 


18 


43 


(+0.08, - 





06) 






13 


59 


(+0.05, 


-0 


06) 


14 


02 


(+0.05, - 





06) 


42.1 (+ 11.8) 


78.9 


(±2 


4) 


18 


40 


(-1-0.08, - 





07) 





9742 


13 


78 


(+0.04, 


-0 


05) 


14 


09 


(+0.05, - 





06) 


67.7 (± 6.9) 


75.9 


(±1 


8) 


18 


59 


(-1-0.05, - 





04) 






13 


74 


(+0.04, 


-0 


05) 


14 


03 


(+0.05, - 





05) 


32.3 (+ 10.3) 


75.7 


(±2 


3) 


18 


56 


(+0.05, - 





04) 





9887 


14 


13 


(+0.06, 


-0 


08) 


14 


54 


(+0.07, - 





09) 


71.9 (± 4.9) 


78.3 


(±2 


4) 


18 


95 


(+0.09, - 





07) 






14 


11 


(+0.07, 


-0 


08) 


14 


55 


(+0.07, - 





08) 


46.1 (± 9.4) 


78.9 


(±3 


0) 


18 


93 


(+0.10, - 





09) 





9898 


14 


24 


(+0.09, 


-0 


11) 


14 


51 


(+0.08, - 





10) 


64.1 (± 4.9) 


75.0 


(±3 


2) 


19 


05 


(+0.13, - 





10) 






14 


18 


(+0.08, 


-0 


10) 


14 


54 


(+0.09, - 





11) 


39.3 (+ 11.4) 


77.3 


(±4 


3) 


18 


99 


(+0.13, - 





10) 





0122 


13 


93 


(+0.08, 


-0 


09) 


14 


27 


(+0.08, - 





10) 


51.4 (± 5.0) 


76.6 


(±3 


0) 


18 


74 


(+0.10, - 





07) 






13 


92 


(+0.07, 


-0 


08) 


14 


19 


(+0.07, - 





08) 


26.2 (+ 12.9) 


75.0 


(±3 


6) 


18 


73 


(+0.08, - 





05) 





0147 


13 


92 


(+0.06, 


-0 


07) 


14 


12 


(+0.08, - 





09) 


55.6 (± 5.3) 


73.3 


(±2 


8) 


18 


75 


(+0.11,- 





09) 






13 


87 


(+0.06, 


-0 


07) 


14 


07 


(+0.07, - 





08) 


29.5 (+ 13.6) 


73.4 


(±3 


4) 


18 


69 


(+0.07, - 





05) 





0175 


13 


72 


(+0.09, 


-0 


12) 


13 


95 


(+0.12,- 





16) 


45.2 (+ 12.4) 


74.1 


(±4 


8) 


18 


54 


(+0.19, - 





13) 






13 


67 


(+0.09, 


-0 


11) 


13 


81 


(+0.15,- 





23) 


34.7 (+ 24.1) 


72.0 


(±6 


2) 


18 


50 


(+0.17, - 





12) 





0209 


13 


89 


(+0.04, 


-0 


04) 


14 


28 


(+0.05, - 





05) 


101.2 (+ 6.7) 


78.0 


(±1 


5) 


18 


70 


(+0.05, - 





05) 






13 


89 


(+0.04, 


-0 
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06) 



to be described by the physics of a low density plasma 
with Tc = 78 000 K we obtain ji = 0.64 [for solar com- 
position according to Anders & Grevesse (1989)]. If we 
identify the stochastic velocity of 68 kms~^ with the tur- 
bulent velocity we obtain a scale height of 3.4 x 10^^ cm, in 
agreement with the best-fit parameter of the hydrostatic 
model. However, at the impact parameter of p = 5.2 Rg 
the theoretical column density exceeds the observed up- 
per hmit by a factor of 3 (see Fig. ||). As expected the 
hydrostatic equation provides an adequate description of 
the inner shell {p < 3 Rg), but overestimates the densities 
in the outer envelope considerably. 

3.3.2. Kinematic model: /3 power-law 

A more realistic model of the large-scale density distri- 
bution requires an appropriate kinematic approach. We 
follow two different strategies to introduce a velocity law: 
an analytical (3 power-law approach traditionally used in 
C Aurigae studies, and a hydrodynamic description of an 



isothermal, isoturbulent wind (where P/p = const). For 
the sake of simplicity we adopt a steady, spherically sym- 
metric outflow. We are aware of the fact that these as- 
sumptions may be a point of criticism, since the obser- 
vations suggest a considerable variability. However, the 
amplitude of these variations are small and justify a kind 
of " average" description. 

The first kinematic model uses a velocity law proven to 
be adequate for different classes of stellar winds including 
the envelopes of evolved late-type stars. The analyses of 
several Q Aur binaries suggest a velocity field of the form 

v{r)=voo{l-^y , (7) 

where Voo is the terminal wind velocity and (3 the accel- 
eration parameter which defines the steepness of the spa- 
tial velocity gradient associated with the expansion. The 
wind velocity v(r) increases monotonically outward and 
approaches the terminal velocity at large distances from 
the giant. The number density is assumed to vary in accor- 
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dance to the equation of continuity. The hydrogen density 
becomes 



nH(r) = np 



RAW r 



(8) 



where the density parameter rt^ is related to the mass-loss 
rate and terminal wind velocity: 



M X 



H 



(9) 



The hydrogen abundance by weight X-^ is fixed at 0.70 as- 
suming solar chemical composition. The least square pro- 
cedure yields a set of best-fit parameters np ~ 6.6 ± 0.7 x 
10^ cm~3 and 13 = 0.60 ± 0.14. 

The small value of (3 indicates a rapid acceleration to 
~ Q.bVoa within 1.5 Rg. We note that the wind mod- 
els of "classical" Q Aurigae systems (C Aur, 31 Cyg , 32 
Cyg) indicate a gradual ac celeration with (3 in the range 
between 2.5 and 3.5 (e.g., Schroder 1985 ). In the case of 
HR 6902 the terminal velocity eludes a spectroscopic de- 
termination, since we can only observe the tangential ve- 
locity component. The poor quality of the spectra does 
not allow the use of alternative velocity indicators like 
Mgii absorption. As a consequence the mass- loss infor- 
mation is restricted to the ratio — . The winds from hy- 
brid stars show high terminal velocities (typically 50 to 
150 knis^^) as measured by the circu mstellar Mgii ab- 
sorption (e.g., Duprec fc Reimers 1987 ). Furthermore, the 
flow velocities of luminous cool stars are always less than 
the escape speed from the photosphere (for HR 6902 we 
find Vcsc = 211 kms^^). So we expect Voc to lie in the 
range 50 — 200 kms^^. The corresponding mass-loss rate 
is 8 X 10^^^ — 3.4 X 10^^^ Mq yr~^ with a statistical error 
of - 10%. 



3.3.3. Kinematic model: P p wind 

Finally, we try to attack the mass motion using the equa- 
tions of gas dynamics. Conservation of momentum gives 
the one-dimensional Euler equation 



dv 
dr 



p dr 



(10) 



where /(r) denotes a potential extra force per unit-mass. 
In the present study we set f — and introduce an ef- 
fective pressure term containing an unspecified nonther- 
mal contribution assuming a proportionality of the form 
Peff oc p. This simplification leads to a Parker-type wind. 
The Euler equation can be integrated immediately and 
yields a solution with a positive velocity gradient at all 
distances (e.g., Mihalas 1978): 



■ exp 



exp 



2rc 
r 



(11) 



19.0 



18.5 



o 




18.0 - 



17.5 



Fig. 5. Comparison of best-fit density models with empir- 
ical hydrogen column densities (dotted: hydrostatic ap- 
proach; solid: P power-law; dashed: P (x p model) 

The transcendental equation can be solved numerically to 
give v{r)/vc- The critical radius is defined by the singu- 
larity of the Euler equation: 



GMg 
'2^ 



(12) 



The hydrogen density may be expressed by an equation of 
the form 



nu{r) = Tic 



Rs 



v{r) ' 



(13) 



where the density parameter ric is related to the mass-loss 
rate by 



M X 



H 



(14) 



This formalism suggests varying the density parameter 
and the critical velocity to model nY{_{r). The critical ve- 
locity specifies the singularity of the momentum equation 
and thus determines the slope of the resulting density 
distribution. We note that there is no need to take ex- 
plicitly care of specific pressure constituents. Our least 
square fit procedure yields Uc — 1.5 ± 0.3 x 10^ cm^'^ and 
Vc = 89.5 ± 13.7 kms~^. The density parameter leads to 
a mass-loss rate of M = 3.4 ± 1.2 x 10"" Mq yr-^ 

Using Eq. |l2| the critical radius becomes ~ 1.4 Rg, 
indicating a high velocity close to the surface of the gi- 
ant. Fig. ^ demonstrates that both wind models lead to 
similar density laws. Indeed, the approximate relationship 
between Vc and [3 derived by Harper et al. [1995, Eq. (10)] 
yields an acceleration parameter of 0.54 quite similar to 



8 



T. Kirsch et al.: The warm circumstellar envelope and wind of the G9 lib star HR 6902 



the result of our beta-law analysis. This finding may jus- 
tify the P power-law approach in other red giant winds, 
where we have an unknown hydrodynamic situation. 

Assuming the effective pressure to be supplied by both 
gas and turbulence pressure the critical velocity is given 

by 



/ kT 



(15) 



The best-fit model with Vc ~ 89.5 kms^^ yields a turbu- 
lence of Vturh = 119 kms~^ which exceeds the observed 
stochastic velocity by a factor of 1.8. This discrepancy 
may become even worse when we allow for macroscopic 
line broadening. Either the identification of turbulent and 
stochastic velocity is completely inappropriate, or the mo- 
mentum Eq. ^ contains unknown external forces or pres- 
sure terms. Furthermore, the relation between the ob- 
served line broadening velocity and the turbulent velocity 
is not clear and depends on the driving mechanism(s), the 
geometry, and the statistics of the absorbing structures. 

In the preceding analysis we have ignored line broad- 
ening due to the wind velocity. In order to examine the 
effect on the derived parameters we repeat the absorption 
line analysis with the complete optical depth quadrature 
according to Eq. |2[ The wind outflow is assumed to be 
well described by the /? power-law. Again we consider ter- 
minal velocities ranging from 50 to 200 kms^^. It turns 
out that high velocity models with Voo ^ 150 kms~^ can- 
not reproduce the observed profile shapes adequately. The 
deduced stochastic component of the line broadening lies 
in the interval 20 — 100 kms^^ showing partly a large 
standard error. The range of stochastic velocities reflects 
both noisy lUE spectra and true irregularities of the local 
velocity field. As an example we present in Table |^ the 
results for Voo = 100 kms~^ which represents the most 
probable fit model. Though this model yields consider- 
ably reduced stochastic velocities with a weighted mean 
(i.e., averaged over all heights) of 47 ± 16 kms~^, the de- 
rived column densities are fairly close to the quasi-static 
model. This finding even holds for larger wind speeds and 
we conclude that the quasi-static analysis has produced 
column densities with a reasonable degree of confidence. 

4. Emission line diagnostics 

Reimers et al. (1990a) have shown that the CiV and Siiv 
emission lines are formed by resonance scattering of B star 
photons in the envelope of the primary. Profile modeling 
requires the solution of a two-dimensional radiative trans- 
fer problem. The difficulty arises due to the fact that the 
light source (i.e. the B star) is displaced from the center of 
wind symmetry. It is well-known that the scattering lines 
are formed in a large volume compared to the scale of the 
binary orbit ( Baade et al. 199^ ). Hence, the emission line 
diagnostics allows to study the large-scale properties of 
the circumstellar envelope. 



For the 2-D radiative transfer calculations we employ 
the SEI (Sobolev with exact integration) method, i.e., the 
source function is approximated using the escape prob- 
ability formalism followed by an exact formal solution 
using short characteristics. This procedure is simple in 
concept and can be applied to very extended geometrical 
grids. In normal C Aur systems the Sobolev approach leads 
to erroneous results, since the wind velocities are typi- 
cally only a few times the stochastic velocities. However, 
a com parison with a more refined 2-D radiative transfer 
code ( Baade 1989 1990) demonstrates the reliability of 



the Sobolev solution in the case of HR 6902. The reason is 
the low mass-loss rate (M w 10~^^ M0 yr~^) compared to 
the "classical" C Aur systems (M w IQ-S-IQ-** Mq yr^^). 
As a consequence the line opacities are small and multiple 
line photon scattering is of minor importance. 



4-1- Sobolev approximation 

In the framework of the first order escape probability the- 
ory the scattering integral is given by (e.g., Mihalas 1978[ ) 



(J) = [l-/3(r)]^(r)+/3e(r)/c 



(16) 



where Ic denotes a constant core intensity representing 
the B star. The escape probabilities are given by 



Pir)-- 

and 

/3c(r) 



47r J 4^ 



exp ( 



■dfi, 



47r 



1 - exp (-Ts) 



(17) 



(18) 



with the Sobolev optical thickness Ts = kon{r)L{r, n). The 
so-called Sobolev length in the direction n may be defined 
in tensor notation: 



L(r,n) 



Edui 
ntTij — 
or, 



(19) 



The quantity /3(r) denotes the probability that a photon 
emitted at the point r escapes from the outer boundary 
and /3c (r) symbolizes the probability that a photon emit- 
ted by the core (i.e. the B star) penetrates to the point 
r. In the case of pure resonance scattering the line source 
function is simply given by 



S{v) 



I3cir)lc 



(20) 



Detailed information can be found in Hempe (1982) who 
applied the Sobolev formalism to the specific conditions 
of the binary configuration. With known source function 
on a two-dimensional grid, we can calculate the emergent 
radiation field. The intensities are supplied by the formal 
solution of the transfer equation and can be integrated 
along a set of parallel rays through the envelope. In a 
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Table 4. Observed emission line fluxes 



Ion 


A 


f 




(A) 


(fO"^'' erg cm-^s-i) 


Si IV 


f395/f408 


14.2 ±3 


Civ 


f548/f550 


24.7 ±4 



subsequent step the integration over all solid angles and 
normalization yields the residual flux F^/Fc. For a detailed 
description we refer to Baade (1990) and Kirsch & Baade 
(1994). 

4.2. Application to HR 6902 

To calculate synthetic emission lines we extrapolate the 
wind models derived in Sect. ^ to the whole envelope as- 
suming the temperature and thus the ionization state to 
remain fixed. Though this simplistic approach implies spu- 
rious energy sources to maintain the temperature it may 
be sufficient to explain qualitatively the observed flux. 

In the first step we use the C iv and Si iv line ratio ob- 
served during total eclipse to derive the ionization temper- 
ature. Unfortunately, the clearly visible C ii line eludes the 
diagnostic procedure, since the photospheric background 
profile is not known precisely. Furthermore, the loci of line 
formation may be quite different. We calculate a grid of 
theoretical flux ratios in the temperature range between 
60 000 and 90 000 K. The ionization problem is solved us- 



ing the CLOUDY code as described in Sect. 3.2. To con- 
sider the influence of different wind velocities we have a 
closer look at four different models defined either by the 



(3 power-law or the P p wind according to Sect. 3.3. 

As shown in Fig. ^ the theoretical flux ratios de- 
pend only slightly on the specific wind model. The ob- 
served line ratio suggests an electron temperature of = 
82 500 ± 3000 K consistent with the mean temperature ob- 
tained on the basis of the absorption line analysis. We use 
the derived electron temperature to calculate integrated 
line fluxes of the C iv and Si iv emission lines. The ob- 
served values (Table ^) can be reconstructed using the 
parameters np = 2.9 x 10^ cm~'^ for the (3 power-law and 
ric = 8.3 X 10^ cm~^ for the P oc p wind model. In both 
cases the density parameter differs only by a factor of ^ 2 
from the results of the absorption line analysis. The dis- 
crepancy may be a result of the inadmissible extrapolation 
of the isothermal density distribution. Indeed, we would 
expect a decreasing temperature at large distances from 
the giant causing a reduced absorber density. 

The presence of a macroscopic velocity field will in- 
fluence the profile shapes and cause wavelength shifts of 
the emission lines. However, the low resolution lUE spec- 
tra (i? w 6 A) does not allow to examine these effects 
more closely. To demonstrate the expected situation we 
present in Fig. a selection of theoretical flux profiles. It 



> 

> 
O 



max 




observed ratio 
min 



T (10^ K) 



Fig. 6. Civ/Siiv line flux ratio vs. electron temperature. 
The observed ratio is obtained from the total eclipse spec- 
trum SWP 39795. The theoretical curves are presented for 
the P power-law approach with different parameter sets 
[M(MQyr-i),Doo(kms-i)]: Solid [8 x 10-12,50], dashed 
[1.7 X 10-", 100], dotted [3.4 x 10-",200]. The P oc p 
model is indicated by a dashcd-dotted curve 




1548 1549 1550 1551 
Wavelength (A) 



1552 



Fig. 7. Theoretical line fluxes of the C IV resonance dou- 
blet during total eclipse {(p — 1.000). The line-center wave- 
lengths are indicated by vertical lines. We present emer- 
gent profiles for the same wind models as shown in Fig. |^ 

is obvious that the asymptotic flow velocity determines 
the character of the line profiles. In combination with a 
semi-empirical density model the line positions and pro- 
file shapes give reliable information about the global wind 
structure. This allows us to predict that high-resolution 
UV observations will unravel the kinematic nature of the 
circumstellar envelope of HR 6902. 
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5. Summary and discussion 

In the present paper we have analyzed lUE observations 
of the C Aurigae type binary system HR 6902. The hot 
secondary serves as a convenient probe for the extended 
envelope of the bright giant. By means of a least square 
profile-fit procedure we have deduced CiV and Siiv col- 
umn densities and stochastic velocities. The density distri- 
bution of the inner envelope (< 3 Rg) can be described by 
a quasi hydrostatic model. We have shown that pressure 
support by turbulence is sufficient to explain the empirical 
scale height. However, in order to match the empirical col- 
umn densities adequately we had to introduce an appro- 
priate outflow model. We have demonstrated that both the 
common (3 power-law and a P oc p wind yield appropri- 
ate density distributions and explain the absence of C iv 
and Si IV absorption at larger distances from the giant. 
Fig. H demonstrates that the hydrostatic solution as well 
as the different outflow models lead to similar density laws. 
Despite this similarity the hydrodynamic implications are 
completely different. The required driving force per unit 
mass f{r) should be compared to theoretical models. 

We have derived mass-loss rates in the range 0.8 — 
3.4 X 10~^^ Mq yr"-'^. This result indicates that the mass- 
loss rate of HR 6902 scales roughly with the solar mass 
flux, i.e., M « M0(Rg/RQ)^. Finally, we have demon- 
strated that the emission line spectrum obtained during 
total eclipse confirms essentially the wind models derived 
from the absorption line analysis. This finding supports 
the picture of a steady and spherically symmetric flow 
without excessive deviations from the mean wind model. 

It is shown that HR 6902 possesses warm (To « 
80 000 K) circumstellar material, observed in a range of 
impact parameters between l.I and 3.3 Rg. Recently, 
the innermost region immediately above the photosphere 
{p < 1-1 f^g) has been examined by Schroder, Marshall, 
& Griffin (1996) using optical spectra. They derived a 
temperature of ~ 5000 K and a turbulent velocity of 
^^turb = 15 kms~^. However, in order to fit the Can K-line 
wing they postulated a second absorption component with 
a considerably larger turbulent velocity of ^ 45 kms^^ 
quite similar to the stochastic velocities derived in the 
present work. Combining these results we suspect a dra- 
matic change of the physical conditions in the giant's at- 
mosphere at p « 1.1 Rg associated with a steep rise of the 
temperature from about 5000 to 80 000 K. 

The large nonthermal line broadening parameter and 
its relative constancy (apart from local fluctuations) may 
be a characteristic property of the wind of HR 6902. The 
simplistic identification of turbulence and the observed 
stochastic velocity cannot explain the large pressure re- 
quired to explain the density distribution by a P oc p 
wind. However, the true relation between the observed 
line width and turbulent motions associated with wave- 
driven winds depends on the specific driving mechanism. 



250 




E 10= 
o 



10" 



P power-law 

P cc p 
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Fig. 8. Velocity (a), density structure (b), and nonthermal 
momentum deposition (c) plotted as a function of radial 
distance. The (3 power-law models are labelled with the 
corresponding terminal velocity 

geometrical details, and the randomness of the absorbing 
structures (e.g. 



ludge 1992) 



The wind models suggest a rapid wind acceleration 
within ~ 0.5 Rg above the photosphere. This result is sim- 
ilar to the outflow characteristics of hybrid stars [a Aqr 
(G2 lb), 13 Aqr (GO lb), 9 Her (Kl Ha), l Aur (K3 II) 



7 Aql (K 311), and g TrA (K3 II); [Brosius fc MuUan 1986 
pupree ct al. 199^ ; Harper et al. 1995 1. Furthermore, a 



recent analysis of ORFEUS II spectra of the hybrid stars 



a TrA and a Aqr (Dupree & Brickhouse 1998) have re- 



vealed that the stellar winds are not confined to cool cir- 
cumstellar material. They pointed out that ions like Sim 
through O VI indicate expanding material up to a tem- 
perature of 3 X 10^ K with wind velocities from 100 to 
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200 kms ^. Wc argue that the wind properties of hy- 
brid stars are comparable to the extended atmosphere of 
HR 6902. 

With its location in the HR diagram HR 6902 
represents the physically important connection between 
solar-like stars with hot coronae and late type supergiants 
with cool outer atmospheres. The binary technique 
provides spatially resolved information about the phys- 
ical structure of the extended envelope and wind that 
arc unobtainable by other means. With lUE we have 
gained a first insight, but additional high resolution UV 
observations are indispensable to construct a detailed 
model. With this unique opportunity HR 6902 may play 
a crucial role in the elucidation of atmospheric heat- 
ing processes and mechanisms to drive winds of cool stars. 
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